Surface polar optical phonon interaction induced many-body effects and hot-electron 

relaxation in graphene 
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We theoretically study various aspects of the electron-surface optical phonon interaction effects 
in graphene on a substrate made of polar materials. We calculate the electron self-energy in the 
presence of the surface phonon-mediated electron-electron interaction focusing on how the linear 
chiral graphene dispersion is renormalized by the surface phonons. The electron self-energy as well 
as the quasiparticle spectral function in graphene are calculated, taking into account electron-polar 
optical phonon interaction by using a many body perturbative formalism. The scattering rate of free 
electrons due to polar interaction with surface optical phonons in a dielectric substrate is calculated 
as a function of the electron energy, temperatures, and carrier density. Effects of screening on the 
self-energy and scattering rate are discussed. Our theory provides a comprehensive quantitative 
(and qualitative) picture for surface phonon interaction induced many-body effects and hot electron 
relaxation in Dirac materials. 



I. INTRODUCTION 

Lattice vibrations are inevitable intrinsic sources of 
carrier scattering in semiconductors and often domi- 
nate transport near room temperatureii2,. In general, 
scattering by acoustic phonons is the main contribu- 
tion to the electronic resistivity at low temperatures, 
and the longitudinal-optical (LO) phonons in polar ma- 
terials limit the room temperature mobility through 
their coupling with carriers via the long-range Frohlich 
interaction^. Since graphene is made of carbons it is 
a non-polar material and the lack of strong long-range 
polar optical phonon scattering leads to very high in- 
trinsic (i.e., phonon-limited) room-temperature graphene 
mobility^""— and indeed the room-temperature graphene 
mobility is typically very high in disorder-free suspended 
graphene sample a^^d^ . The weak deformation poten- 
tial scattering from the thermal lattice acoustic phonons 
mainly limits the intrinsic mobility of graphene at room 
temperature and the acoustic phonon scattering gives a 
quantitatively small contribution even at room temper- 
ature due to the high Fermi temperature at high carrier 
density (and the small deformation potential coupling) 
of graphene^"—. 

However, most current available graphene samples for 
fundamental studies and certainly for technological appli- 
cations use a polar substrate such as Si02, SiC, or Hf02. 
Such polar substrates allow for the existence of polar op- 
tical phonons localized near the graphene-substrate in- 
terface which could be an important scattering source 
for graphene carriers through the long-range Frohlich 
coupling. The general eigenvalue equation for a surface 
phonon mode has been derived at a planar interface be- 
tween two polar dielectrics^"^"^^. Since the surface op- 
tical (SO) phonon is a well-defined surface property at 
the interface of the polar semiconductors^^' it is pos- 
sible that carriers in graphene layer couple to the SO 
phonon of the underlying substrate lattice via the long- 



range polar Frohlich coupling. Recently, the possible role 
of surface optical phonon on high temperature graphene 
mobility has been pointed out^"^. Due to the polar na- 
ture of the gate dielectrics used in graphene field effect 
transistors, the carriers in the conducting channel elec- 
trically couple to the long-range surface phonons created 
at the dielectric interface. In general, the surface optical 
phonon contributes little to the carrier resistivity at low 
temperatures due to its high energy {losop ^ bOmeV). 
However, the contributions of other scattering processes 
to the graphene resistivity are small at room temper- 
ature and therefore the SO phonon may become the 
dominant scattering mechanism in graphene on a polar 
dielectric substrate at high temperatures. It is there- 
fore possible, perhaps even likely, that in the context 
of graphene technological applications, SO phonon ef- 
fects are particularly important and must be taken into 
account. In addition, in a recent experimenti2r— , the 
strongly coupled plasmon-phonon mode dispersion has 
been measured by the angle-resolved reflection electron- 
energy-loss spectroscopy and it is found that the electron- 
phonon coupling give rise to the many interesting many- 
body effects in graphene on a polar dielectric substrate^i. 
This again points to the possible importance of carrier- 
SO phonon coupling in determining the electronic prop- 
erties of graphene on a substrate. 

Since the electron-optical phonon interaction leads to 
many-body renormalization of the single-particle free 
carrier properties^i^ while at the same affecting the 
transport properties^ a good understanding of electron- 
SO phonon coupling is thus important in developing 
quantitative theories for many different experimental 
studies in graphene. In this paper we calculate the elec- 
tron self-energy and spectral function of graphene on a 
polar substrate such as Si02, SiC, or Hf02. Our goal in 
this paper is to provide a fairly detailed theoretical pic- 
ture for the electron-SO phonon interaction induced polar 
scattering of free carriers. Many electronic properties of 



2 



graphene on polar substrates are expected to be modified 
by electron-SO phonon scattering via the Flohlich inter- 
action. On a microscopic level this interaction renormal- 
ize the Fermi velocity of graphene. A thorough quantita- 
tive understanding of the electron-SO phonon interaction 
induced many-body effects in high mobility graphene is 
very important with respect to both the device opera- 
tion and the fundamental physics of graphene on sub- 
strates. This surface polaronic (i.e. interaction with SO 
phonons) effect is an important factor in determining the 
high temperature mobility and the hot electron proper- 
ties of graphene on substrates. 

In this paper, we focus on the effect of electron-SO 
phonon interaction effects on the real, Re[I]], and the 
imaginary part, Im[E], of electron self-energy. From 
Im[I]], we can extract the quasiparticle scattering rate 
(scattering lifetime), which gives information relevant to 
possible applications, such as the hot electron relaxation 
rate and the energy dissipation rate of injected carriers 
in a graphene-based device. We calculate the scatter- 
ing rate as a functions of wave vector, energy, density 
and temperature and discuss the role of screening ef- 
fects. Since the electron- surface polar optical phonon 
interaction has a long range nature the consideration of 
screening effects is crucial in order to understand quan- 
titative physical properties of graphene. We find that 
the scattering rates with screened interaction is reduced 
by a factor of five, compared with the results with un- 
screened interaction. Our consideration of screening (i.e., 
the static random phase approximation (RPA) screen- 
ing) may overestimate the screening effect on scattering 
rate at low temperatures, but it provides an upper limit 
of the screening effects. In addition to calculating the 
properties of Im[I]] due to SO phonon interaction, we 
also calculate the Re[I]], which gives us the correspond- 
ing graphene velocity renormalization due to the surface 
polaronic effect. We also provide results for the graphene 
spectral function due to electron-SO phonon scattering, 
which can be directly observed in scatting tunneling spec- 
troscopy (STS) and angle-resolved photo-electron spec- 
troscopy (ARPES). Our theory and results are directly 
relevant to the electronic properties of doped or gated 
graphene on polar substrates where SO polaron effects 
may very well be significant. 

The paper is organized as follows. In Sec. II the gen- 
eralized theory is presented to calculate the self-energy 
of electron in the presence of electron-SO phonon inter- 
action. Sec. Ill presents the results of self-energy and 
spectral function. We also show the scattering rate as 
a function of energy, momentum, density, and tempera- 
ture. We summarize in Sec. IV with a discussion. 



II. THEORY 

The Hamiltonian of graphene is well- approximated by 
a 2D Dirac equation for massless particles, Hq = vpcr ■ k, 
where vp is the Fermi velocity of graphene, a is the pseo- 



dospin 2x2 Pauli matrices, and k is the 2D momentum 
relative to the Dirac points (we use h = 1 throughout this 
paper). The two components of the spinors correspond to 
occupancy of the two sublattices of the honeycomb struc- 
ture in a hexagonal lattice. This Hq gives a linear energy 
dispersion relation Sk.s = 5WF|fe|, where s = +1 (— 1) 
for the conduction (valence) band. The corresponding 
density of states (DOS) is D{e) = gsgvl^l/ {2nvp), where 
gs = 2, gy = 2 are the spin and valley degeneracies, re- 
spectively. 

The central quantity we calculate in this paper is 
the leading order electronic self-energy correction due 
to electron-surface optical (SO) phonon interaction via 
Frohlich polar couplingi^. In the presence of the long 
range electron-SO phonon coupling, electrons interact 
among themselves through the Coulomb interaction and 
through virtual- SO phonon exchange via the Frohlich in- 
teraction. The electron-SO phonon interaction is given 

by 



ph 



kq ss' 



where c^^ is the electron (s 



-1) or hole (s 



(1) 



-1) 



creation operator, 6j and bq are creation and destruction 
operators of surface optical phonon, and the interaction 
matrix element MJ"* is defined by 



where F^k is the chiral spinor and given by 



1 



(2) 



(3) 



with s = ±1, 0k = tau.~^ (ky/kx) is the angle between 
fc, k' , arises from the overlap of \sk) and |s'fe').-2. The 
strength of the coupling M (q) is given by 



where 



[Miq)? 



7re ujso 



1 eo + l 



(4) 



(5) 



and e(g) is the dielectric (screening) function of 
graphene^, d is the separation distance between 
graphene layer and substrate, ujso is the surface opti- 
cal phonon frequency, and eo (coo ) is the static (high fre- 
quency) dielectric constant. In Eq. (jH), we set e(g) = 1 for 
the unscreened case. For screened interaction we use the 
static RPA dielectric function appropriate for graphene, 
which is give by 



eiq) = l + viq)n{q), 



(6) 



where v(q) = 2'Ke^ j nq is the electron-electron Coulomb 
potential with the effective background lattice dielec- 
tric constant k = (eo -|- l)/2, and ^{q) = \\{iq,T) is 



3 



the electronic polarizability function of graphene which 
depends on both wave vector and temperature. For 
isotropic media the frequency of SO phonons to so is re- 
lated to the transverse optical (TO) bulk phonon loto els 
Loso/uJTO = A/(eo + l)/(eoo + 1)."'^^ Note that the bulk 
longitudinal optical phonons uj^o a-nd ujto are connected 
with the dielectric const ants by the Lyddane-Sachs- Teller 
relation^ wlo/wto = i/eo/eoo- 

The SO-phonon mediated electron-electron effective 
interaction (i.e. in the leading-order screened Frohlich 
coupling) is given by 



vtl(q,u)^[AUq)]^D^{u:)F^l,, 



(7) 



which depends on both wave vector and frequency. In 
Eq. ([7|) Do{uj) is the unperturbed SO-phonon propagator 

I 



and given by 



and 



TPSS 

^kk' 



Doicu) 



^Ik-^s'k+q 



(8) 



""SO 



[l + ss'cos{6kk')]/2. (9) 



The self-energy E within the leading-order approxima- 
tion is give by^ 

s' q,'ii'n 

(10) 

Here, /3 — l/fcsT, s, s' — ±1 are band indices, 
Go,s(k, jaj„) — l/(jw„ — Cks) is the bare Green's function, 
LOmVn are Matsubara fermion frequencies, and Af^^ the 
interaction matrix element. After the standard analytical 
continuation from ia;„ to w -I- iO"*" the retarded self-energy 
is obtained, 



Sf (fe,c.) 



E 



-2qd 



dq M^e 
(27r)2 qe(q) 



-F,s>{k,k~q) 



No + npiks') 



No + 1 - npiks') 



oj + ir] + Luso - £.k's' LO + irj-LOso-ik'. 



(11) 



where s = VF\k\ — /i, is the electron energy relative to 
the chemical potential ^ {— Ep at T = 0). The func- 
tion A^o and np are Bose and Fermi distribution function 
defined by 



and 



np{ks) = 



1' 



(12) 



(13) 



r 



where the chemical potential ^ at finite temperature is 
determined from the total electron density n of the sys- 
tem. The self energy is an explicit function of the in- 
dependent variable momentum (fc) and energy u. We 
obtain the real and imaginary parts of the electronic self 
energy 



Re[I],(fe,^)] ^Mo^X J l^^-'^'^'^') 



NQ + np{^k's') , No + l- np{S,k's') 



i^so - i 



k' s' 



^ — ^SO — £,k's' 



(14) 



/k' dk' 
^^C-(fc,fc')[(^^0+"F(&'.'))^(^ + ^SO-&'.') + (^0 + l-nF(&'s'))'5(t^-t^SO-a's')], 

(15) 



where 



Gss'{k,k') = j de- 



-■2\k-k'\d 



k - k'\e{\k - k'\) 



—Fss'{k,k'), (16) Eq. p^. becomes 



I 

where 6 is the angle between k and k' . 

For T = the imaginary part of the self energy. 
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{lu + u>so + Ef)]s'{lo + ujso + Ep)9{—Lo — u)q)9[s' (lo + + -Ef)] 

s'=±l 

+Gss' [fc, s'(w - wso + Ef)]s'{lo - ojso + Ef)0{uj - uo)0[s'{u - ojso + Ep)]} , (17) 



where A = -ttM§ /{2'k)'^vf and 9{x) = 1 for x > and 
for X < 0. The first term in Eq. ^ denotes SO 
phonon absorption by the electron, whereas the second 
term denotes the emission of SO phonon. From Eq. (jl7p 
we find that the phonon emission is only ahowed for 
Lu > uJso, and the phonon absorption for ui < —loso- 
Thus, due to the phase space restriction there is a re- 
gion where lm\^{kF ^oj)] becomes zero for energies within 
UJSO of the Fermi energy, i.e., for —ivso < < uJso 
Im[I](fc_F,a;)] = 0. Note that uj — corresponds to the 
Fermi energy. The vanishing of Im[I]] indicates that a 
quasiparticle within this region cannot decay by emitting 
(or absorbing) surface phonons because there is no avail- 
able final state to decay. A quasiparticle with an initial 
state energy w has a final state energy uj — ujso by emit- 
ting an SO phonon. The final state energy of the electron 
must greater than the Fermi energy (i.e., lu—ujso > 0) for 
a decay process to occur. Otherwise the transition can- 
not occur because the final states are all occupied. Thus, 
for the quasiparticle with energy lo < ujso, Im[E] = 0. 
An electron outside this energy range can lose energy by 
emitting SO phonons. For a quasiparticle with the initial 
energy uj has the final state lo -\-ujso by absorbing an SO 
phonon. Thus, the quasiparticle with the initial energy 
—UJso < ^ cannot decay by absorbing an SO phonon. 
Thus, for —UJso < uJ < uJso, Im[I](feF,w)] = 0. In addi- 
tion to this forbidden region, the decay of an electron to 
the Dirac point is not allowed because the DOS is zero 
at that point, which occurs at w = —Ef — uJso (i-e., the 
quasiparticle with the initial energy —Ef — ujso cannot 
decay to the final state energy —Ef (i.e. Dirac point) af- 
ter absorbing an SO phonon). However, we note that this 
forbidden gap exists only for the self-energy contribution 
from surface optical phonons and at zero temperature. 
The quasiparticle will have finite values of Im[I]] in this 
energy region arising from scattering by impurities and 
from electron-electron interactions. 

From the self energy, E(A;, w) , we can obtain the single- 
particle spectral function 



A{k,uj) = -2ImG(fc, 



(18) 



where 



[G(fc,w)]-i = [Go(fc,c^)]-i - I](A;,w) (19) 
Thus, we have 

21mT,(k, uj) 

A{k, ij) = J. 20 

[uj-^k-ReJ:{k,u;)]^ + [lmJ:{k,uj)f 

The spectral function A(k,uj) can roughly be interpreted 
as the probability density of the different energy eigen- 
states required to make up a specific k state. The spectral 



function shows a sharp peak at the quasiparticle state 
and becomes a d function when Im[I]] = 0. The spectral 
function must satisfy the sum rule 



—A{k,u) 
Ztt 



(21) 



and in all our numerical calculations we have explicitly 
checked that this sum rule is generally satisfied to within 
less than a percent. 

The inverse quasiparticle lifetime (or, equivalently, the 
scattering rate) rs(fc) of state \sk) is obtained by setting 
the frequency in imaginary part of the self-energy to the 
on-shell (bare quasiparticle) energy £^sk = svF\k\ — Ef, 
i.e. 



r,(fc) = 2Im[I]f (fc,a.)]. 



(22) 



The self-energy approximation used here is equivalent 
to the generalized Born approximation for the scatter- 
ing rate. Note that the integrand of Eq. (ITSI) is non-zero 
only when ^ks = ii^so + Cfc's', which correspond to emis- 
sion (-I-) and absorption (— ) of surface optical phonon, 
respectively. 



III. RESULTS 

In this section we present the results for graphene on 
SiC. The parameters used in this calculation are follow- 
ing: the surface optical phonon frequency ujso = 117.6 
meV, the longitudinal optical phonon frequency, ujlo = 
120.7 meV, the static (high frequency) dielectric constant 
£0 = 9.72 (ecx, = 6.52). 

In Fig. [T] we show the calculated self-energy as a func- 
tion of energy {uj/Ep) for k = kp, d = 0, and T = 0. 
Note that uj = Q and -Ep represent the Fermi energy and 
Dirac point, respectively. The density n = 2 x 10^^ cm~^ 
is used in the calculation, which corresponds to Ep — 165 
meV. The Fermi momentum {kp) and the Fermi energy 
{Ep, relative to the Dirac point energy) of graphene are 
given by kp — {^-ktl/ gsg^)^/"^ and \Ep\ = vpkp where n 
is the 2D carrier (electron or hole) density. (S-) rep- 
resents the self energy for s = 1 (s = — 1). In Fig. [1] (a) 
and (b) the self-energies are calculated without screen- 
ing effect by putting e{q) = 1 in Eq. (fTTj) . but in (c) and 
(d) the RPA screening function is used to calculate the 
self-energies. As discussed in Sec. II, Im[I]s] vanishes 
for —UJso < UJ < UJso and at a; = —Ep — ujso, and the 
range of the forbidden gap is independent of screening. 
For s = — 1 there is a sharp strong peak in Im[I]_]. This 
singular behavior for unscreened self energy for s = — 1 
is unphysical and arises from the long range nature of 



5 





2 

co/Ef 



FIG. 1. Calculated self energies as a function of energy, u/ Er- 
in (a) and (b) E±(fc,aj) is calculated without screening effect 
[i.e. e(q) = 1 in Eq. ([Tl|l and for fe = fe^, d = 0, and T = 0. In 
(c) and (d) the self-energy E±(fc,a;) is calculated with RPA 
screening function. Solid red (dashed black) lines indicate 
imaginary (real) parts of the self energy. 



FIG. 2. Imaginary part of the self energy as a function of 
energy for k — kp and T = and for different values of 
distance between graphene and the dielectric substrate, d = 
(solid line), Inm (dashed line), 2nm (dot-dashed line) (a),(b) 
are calculated without screening effect and (c),(d) with RPA 
screening function. 



the electron- SO phonon interaction, i.e., Vph cx 1/q. The 
RPA screening makes two significant modifications to the 
self energy compared with the results with the unscreened 
electron-SO phonon interaction. The screening effects re- 
duce the magnitude of both Re[i;] and Im[I]] by roughly 
a factor of 10, and removes the logarithmic singularity 
in ImS]_ occurring at a; = —Ep — LOgo ~ vpk. Since 
screening makes the long range electron-phonon interac- 
tion finite, the singularity at q = (or, k' — k) is removed 
by changing q to ge(g) in Eq. p7)) . 

In Fig. [5] Im[S±(fc, cj)] is shown as a function of en- 
ergy Lo for k = kp and T = 0, and for different values 
of d = 0, 1, and 2 nm, where d is the distance between 
graphene and the interface of the dielectric substrate. In 
Fig. (Ha) (b) the self-energy is calculated without screen- 
ing effects and in (c) (d) the RPA screening in electron- 
phonon interaction is considered. As the distance d in- 
creases the magnitude of the Im[S] decreases due to the 
exponential dependence of the interaction, i.e., e"^*'' in 
the electron-SO phonon interaction, [see Eqs. ^ and 
([7])]. In Fig. [2ja) and (b) the logarithmically diverg- 
ing peaks at cli = ±i^so and lo — —ikp + i^so arises 
from the singular behavior of the long range electron-SO 
phonon interaction at g = 0. However, as discussed in 
Fig. [Jl this peaks disappear in the presence of screening 
as shown in Fig. [5] (c) and (d). For comparison of these 
results with the imaginary parts of the self-energy cal- 
culated with non-polar in-plane LO phonon (wlo ^ 200 
meV) arising from the lattice vibration of graphene itself 
which has been discussed in Ref. [l^ and [2^, we con- 
sider the self-energy for d = and the unscreened case 
{t{q) = 1). Then, the main difference between these two 



cases appears in Eq. ([T7| . For non-polar LO phonon we 
have Gss' — tt, and for polar SO phonon we have 

GsAkX)=^-^[K{r)~E{r)] 



kk' 
k + k' 
kk' 



k' 



k + k' 



K{r)+E{r) 



(23) 



where K{r) and E{r) are the complete elliptic inte- 
gral of the first and second kinds, respectively, and 
r = 2kk'/{k'^ + fc'^). For the non-polar LO phonon Im[E] 
increases linearly with \uj\ for larger |a;| because the DOS 
increases linearly with For polar SO phonon the 

long range interaction and the linear energy dispersion 
of graphene DOS compensate each other and the can- 
cellation of these two effects makes the self energy sat- 
urate for large \uj\. The emission (or absorption) of SO 
phonon does not depend on the quasiparticle energy for 
large |w| (our numerical results show that the Im[S] sat- 
urates for |a;| > 'iEp). Thus, high energy hot electrons 
decay more effectively by emitting (absorbing) graphene 
LO phonons than low energy electron. However, for SO 
phonons the damping of the electron by emitting (ab- 
sorbing) SO phonons depends weakly on the energy of 
electrons. For finite d, Im[I]] decreases faster than for 
d = as |a;| increases, then the high energy hot elec- 
trons do not decay effectively by emitting (or absorbing) 
remote SO phonons. 

In Figs. [5]— {5] we show the calculated self-energy 
E(/c,cj) and the spectral function A{k,Lu) for different 
values of fc = 0.5 1.0, and 1.5 kp. Figs. |3] and S] show the 
results calculated with unscreened electron-SO phonon 
interaction for s = +1 (conduction band) and s — —1 
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FIG. 3. Calculated self energy E-t- and spectral function A+ 
as a function of energy for different wave vectors, k — 0.5, 
1.0, 1.5 kp and for d — Inm and T = 0. Note that the ener- 
gies are measured in terms of uso instead of Ep. Screening 
effect is not considered, e(g) — 1. The vertical arrows in- 
dicate (S-function peaks at the quasiparticle energy with the 
spectral weights shown next to the arrows. They are located 
at Lj = -0.447, 0.0, and 0.454 Ep for k = 0.5, 1.0, 1.5 kp, 
respectively. We find that from the calculated spectral func- 
tion the renormalized velocity due to the electron-SO phonon 
interaction is reduced, i.e., vp — 0.88vp approximately for 
these results. 



(valence band), respectively. The bare electron (hole) 
energy, E{k) — zLvF\k\, is modified due to the electron- 
SO phonon interaction and the energy dispersion is given 
by the solution to Dyson's equation for a given fc. 



(24) 



where = VF\k\ — Ep and /i = S+(fcF, 0) is the chemi- 
cal potential. A peak in the spectral function corresponds 
to the solution to Dyson's equation and the quasiparticle 
energy for a given k. For s = -f 1 (Fig. [3]) the spectral 
function has a well defined 5-function peak as long as the 
peak is located inside the forbidden region. The vertical 
arrows in Fig.[3]indicate (5- function peaks at the quasipar- 
ticle energy with the spectral weights shown next to the 
arrows. They are located a.i to — —0.447, 0.0, and 0.454 
Ep for k = 0.5, 1.0, 1.5 kp^ respectively, and these quasi- 
particle peaks consume almost 90% of the total spectral 
weight. The remaining part of the spectral weight goes to 
the incoherent background. We find that from the cal- 



FIG. 4. Self energy E and spectral function A for s = — 1 
calculated with unscreened electron-SO phonon interaction. 
S- and A- are shown as a function of energy uj/luso for 
different wave vectors, k — 0.5, 1.0, 1.5kp and for d = 1 
nm and T — 0. The quasiparticle peaks are not coherent for 
s — —1 and the peaks are broadened by emitting or absorbing 
SO phonons. 



culated spectral function the renormalized velocity due 
to the electron-SO phonon interaction is reduced, i.e., 
Vp = 0.88vp approximately for these results. In Fig.|3]we 
show the self-energy S]_(fc,w) and the spectral function 
A-{k,Lu) for s = — 1 and for different values of fc = 0.5 
1.0, and 1.5 kp. For s = —1 the quasiparticle peaks are 
located outside the forbidden region and they decay by 
emitting (or absorbing) SO phonons. As a consequence, 
the spectral function has a peak with finite broadening 
instead of the (5-function peak. The quasiparticle peaks 
for s = —1 (or in the valence band) are not coherent and 
the peaks are broadened by emitting or absorbing SO 
phonons. 

Figs. [S] and [5] show the results calculated with the 
screened electron-SO phonon interaction for s = -f 1 (con- 
duction band) and s = — 1 (valence band), respectively. 
For s = 4-1 we have the S function peak in the spectral 
function corresponding to the solution to Dyson's equa- 
tion and the quasiparticle energy for a given k. However, 
the energy of the quasiparticle (i.e., the location of the 
peak) for a given wave vector is almost identical to that 
of the non-interaction systems. The peaks are located at 
io = -0.49, 0.0, and 0.49 Ep for k = 0.5, 1.0, 1.5 kp, re- 
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FIG. 5. Calculated self energy E-t- and spectral function A+ 
as a function of energy for different wave vectors, k — 0.5, 1.0, 
l.Skp and for d — 1 nm and T = 0. RPA screening is consid- 
ered in the electron-SO phonon interaction. Arrows indicate 
(5-function peak at the quasiparticle energy. The renormalized 
velocity due to the screened electron-SO phonon interaction is 
bareley reduced, giving Vp = 0.98vf approximately for these 
results. 



spectively, and these quasiparticle peaks consume almost 
all (99%) of the total spectral weight. Just less than 
1% of the spectral weight goes to the incoherent back- 
ground. The renormalized velocity due to the screened 
electron-SO phonon interaction is barely reduced, giv- 
ing Vp = 0.98vf approximately for these results. Since 
screening reduces the magnitude of the imaginary part 
of the self energy as shown in Figs. [1] and [2] the spectral 
weight of the incoherent background is diminished in the 
presence of screening, which gives rise to the increase of 
the spectral weight in the quasiparticle peak because of 
the sum rule. In Fig. |6]we show the self-energy E_(fc,cij) 
and the spectral function A — {k, u) calculated with the 
screened electron-SO phonon interaction. For s — —1 
the quasiparticle peaks in the presence of screening are 
also incoherent and the peaks are broadened by emitting 
or absorbing SO phonons. However, the reduction of the 
magnitude of the self-energy due to screening the broad- 
ened quasiparticle peaks become much sharper than the 
peaks without screening. 

In Fig. [7] we show the imaginary part of the on-shell 
self energy, i.e., the self energy at w = = vpk — Ep. 



FIG. 6. Calculated self energy E_ and spectral function A- 
as a function of energy for different wave vectors, k — 0.5, 
1.0, l.bkp and for d = 1 nm and T = 0. RPA screening effect 
is considered in the electron-SO phonon interaction. 



For < k < kp — uJso/vp electrons decay by absorbing 
SO phonons, while for k > kp +ujso/vp electrons decay 
by emitting SO phonons. Thus, Im[Es(fc,^fe)] = for 
kp — ujso/vp < k < kp + uJso/vp. Inside the forbidden 
gap, the electron does not decay by emitting or absorbing 
a SO phonon. For d = and large wave vectors k we have 
Im[E+(fc,^fe)] oc G++{k,k - uJso){k - ujso) « log(fc) and 
Im[E_(/c,a)] « G-+ik,k~ioso){k~ujso) - 2. The 
Im[S-|_(fc, ^fe) increases logarithmically and Im[E_(A:,^fe) 
saturates for large k. Thus, the damping of electrons by 
emitting SO phonon depends weakly on the energy of 
electrons. However, the Im[E_|_(fc, ^fe)] for the non-polar 
LO phonon increases linearly with wave vector and high 
energy electrons decay effectively by emitting more LO 
phonons^^'^^. Screening effects reduce the magnitude of 
the on-shell energy as shown in Fig. and (d). In 
addition to the reduction due to screening the magnitude 
of the imaginary part of the on-shell self-energy decreases 
with increasing d (the distance between graphene and the 
interface of the dielectric substrate) . 

In Fig. [5] the calculated damping rate, r(fc) = 
2Im[E(fc, ^i:)], is shown as a function of wave vector 
for different temperatures. Here, we assume that the 
temperature of electron and phonon is the same. As 
temperature increases the forbidden gap disappears due 
to the thermal smearing of the Fermi surface, which 
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FIG. 7. The imaginary part of the on-shell self energy, 
Im[E(fc, ^fc)], as a function of wave vector for different d = 
(sohd hnes), 1 (dashed lines), 2 nm (dot dashed lines). The 
calculation is done at T = 0. (a) and (b) are the results 
with unscreened electron-phonon interaction for s = +1 and 
s = — 1, respectively, (c) and (d) are results with screening 
interaction for s = +1 and s — —1, respectively. The forbid- 
den gap, where Im[Es (fc, ^fc) =0], is given by kp —ujso/'vf < 
k < kp + UJSO/VF- 



FIG. 8. Calculated damping rate T{k) — 2Im[E(A;, ^fc)] as a 
function of wave vector for different temperatures, T = 0, 0.2, 
0.4, 0.7, 1.0 Tf (from bottom to top), where Tf = -BF/fcs is 
the Fermi temperature. The carrier density n — 2x 10^^cm~^ 
and d — are used in this calculation, (a) and (b) are the 
results with unscreened electron-phonon interaction for s = 
+ 1 and s — —1, respectively, (c) and (d) are results with 
screened interaction for s = -1-1 and s = — 1, respectively. 



loosens the phase space restriction. Thus, SO phonon 
emission (absorption) is allowed for k < kp +ujso/vf 
(k > kp — LUso/vp)- The quantitative effect of screen- 
ing (compared with the unscreened result) is about a 
factor of five reduction of the scattering rate as shown 
in Fig. [5] (c) and (d). In addition to the reduction of 
the rate, screening effects are particularly important in 
r_(A;). In Fig.|S](b) the peak inside the forbidden region 
at T = develops as temperature increases. The peak 
occurs at k — LOsof^vp and arises from the singular na- 
ture of the long range bare electron-phonon interaction. 
The divergence of the interaction at g = disappears in 
the presence of screening effects and the developed peak 
in scattering rate also disappears due to screening effects. 

In Fig.[3]the calculated damping rate at k = kp, r{kp), 
is shown as a function of temperature for different elec- 
tron densities. The results in Fig. [3] are calculated with- 
out screening effects, but we notice that the screening ef- 
fects reduce the scattering rate approximately by a factor 
of five as shown in Fig. [HI The scattering rate increases 
exponentially with temperature due to the increases of 
phonon population, F cx as temperature in- 

creases. Insets in Fig. [5] show the calculated damping 
rate as a function of density at T = 300K. The den- 
sity dependent damping rate at a given temperature is 
non-monotonic and has a local minimum at the density 
corresponding to Ep = ojso- Thus, the damping rate de- 
creases first as the density increases, and has a minimum 
at n = (lolo/vf)'^ /tt. The damping rate increases as 



the density increases further. For s — —1 the damping 
rate shows more complicated density dependence. The 
local minimum occurs at the same density for s = -1-1, 
but there is a sharp peak at low density. The sharp peak 
in r_ (fcj?) arises from the unscreened electron-phonon in- 
teraction, and disappears when we consider the screening 
effects. 

When the scattering rate F is known, the conduc- 
tivity can be calculated by using Boltzmann transport 
equation^ii 

where D{uj) — gsgyUj/ {2Trvp) is the density of states of 
graphene and / is the Fermi distribution function. In 
Fig.Hniwe show the calculated resistivity {p = cr~^) as a 
function of temperature for different d, the distance be- 
tween the interface of the substrate and graphene. From 
the calculated resistivity the mobility is calculated by 
the simple formula /i — 6.24 x 10^/np cm^/Vs where 
n is the carrier density measured in unit of 10^^ cm~^ 
and p in unit of fl. The calculated resistivity limited by 
the electron-SO phonon scattering shows the exponential 
temperature dependence, which is the characteristic be- 
havior of the optical phonon scattering arising from the 
thermal occupation factor in the scattering rate. The 
exponential temperature dependence is not affected by 
the screening effects as shown in Fig. rTOT b) . In insets of 
Fig. [TU] we show the calculated conductivity as a func- 
tion of density for different d. The density dependent 
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FIG. 9. Calculated damping rate at = /cf as a function of 
temperature for different electron densities, which are shown 
next to the lines in the unit of 10^^ cm~^. The results are 
calculated with d = and without screening effects. Insets 
show the calculated damping rate as a function of density at 
T — 300K. The density dependent damping rate have a local 
minimum at the density corresponding to Ef ~ o^so- 



conductivity shows linear behavior, but the slopes are 
different for low and high density regions. The charac- 
teristic electron density Uc separating these two regions is 
given by the density corresponding to the Fermi energy 
Ep u)so- Below the conductivity increases fast 
but above ric the conductivity increases slowly. The ex- 
ponential temperature dependence of the resistivity and 
the linear density dependence of the conductivity have 
recently been observed in the transport measurement^!^. 



IV. CONCLUSION 

We theoretically study various aspects of the electron- 
surface optical phonon interaction effects in graphene on 
a substrate made of polar materials. We provide a rather 
complete set of numerical results of self-energy, spectral 
function, and damping rate both for the unscreened in- 
teraction and for the static RPA screening interaction. 

We find that the scattering rate with screened electron- 



FIG. 10. Calculated resistivity p — a^^ as a function of tem- 
perature for different d, the distance between the dielectric 
interface and graphene. The carrier density n = 2 x lO^'^ 
cm~^ is used in this calculation. In insets conductivity as a 
function of density for different d is shown. In (a) the screen- 
ing effects is not considered, but in (b) the RPA static screen- 
ing is used. The mobility is calculated by the simple formula 
fi — 6.24 X lO^/np cm'^/Vs where h is the carrier density 
(n = n/10^^cm~^) and p in unit of Q. 



SO phonon interaction is reduced by a factor of five 
and the self-energy by a factor of ten, compared with 
the results with unscreened interaction. The static RPA 
screening used in this calculation may overestimate the 
screening effect on scattering rate, but we believe that 
the static RPA screening provides an upper limit of the 
screening effects. That is, our unscreened and statically 
screened results for the surface polar optical phonon scat- 
tering rate provide two extreme bounds for the magni- 
tude of the scattering rate since static screening is an 
overestimation of the actual dynamical screening. We 
mention that the static screening approximation becomes 
very accurate at high carrier density when the typi- 
cal Fermi energy exceeds the SO phonon energy, which 
may already happen at a rather low carrier density 
of 10"'^^cm~^. We therefore believe that our statically 
screened results should be valid for most graphene sam- 
ples used experimentally. In addition to the reduction of 
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the magnitude of the self-energy the screening effects re- 
move the unphysical singular behavior in the self-energy 
for s = — 1 arising from the long range nature of the 
electron-SO phonon interaction. 

We find that the calculated renormalized Fermi veloc- 
ity due to the unscreened electron-SO phonon interaction 
is reduced (with respect to the bare graphene velocity) 
by about 20%, Vp ~ O.Svp, but for the screened inter- 
action it is reduced by only 2%, Vp ^ 0.98vp at a den- 
sity n = 2 X 10^^ cm~^. We note that the renormalized 
Fermi velocity is enhanced about 20% at the same den- 
sity when the electron-electron Coulomb interaction is 
only considered^^. Since our consideration of screening 
in this paper is the RPA static dielectric function we can- 
not state conclusively that the small change of the renor- 
malized velocity in the presence of screened electron- 
phonon interaction arises from the direct cancellation be- 
tween these two interactions. More careful calculation is 
need to confirm this cancellation by considering the dy- 
namical screening effects treating both the electron-SO 
phonon and the electron-electron Coulomb interaction on 
an equal footing. But it appears the graphene velocity 
renormalization due to the electron SO phonon interac- 
tion by itself is a rather small quantitative effect because 



of screening effect. This is certainly true at high car- 
rier densities (> 10^^ cm~^) where our static screening 
approximation should apply well. We also find that the 
scattering rate for the electron-SO phonon interaction de- 
pends weakly on the electron energy due to the cancella- 
tion between the long range nature of the interaction and 
the linear behavior of the density of states of graphene. 
For the in-plane non-polar LO phonon of graphene the 
scattering rate strongly depends on the electron energy 
and high energy electrons are heavily damped by emit- 
ting LO phonons. Since the scattering rate due to the SO 
phonon decreases further due to the physical separation 
between graphene and interface of the dielectric substrate 
the chance of the hot electron damping by emitting SO 
phonons may be very low in the real graphene samples, 
but the high-temperature transport in clean graphene 
may very well be dominated by SO phonon scattering 
from the substrate. 
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